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High- Order Methods 

Discontinuous Galerkin (DG) methods by Reed and Hill 1973, 
Cockbum and Shu 1990’s, Bassi and Rebay 1997, 2000 ... 

• Integral form, stable, powerful machinery 

• Not intuitive 

Staggered-Grid methods by Kopriva and Kolias 1996; Spectral 
Difference (SD) scheme by Liu, Vinokur, and Wang 2004, ... 

• Differential form, simple and intuitive 

• Mildly unstable 

Flux Reconstruction methods (FR, Huynh 2007, Wang and Gao 
2009, Jameson 2010, Vincent, Castonguay, Jameson 2011, ...) 

• Differential form, recovers DG, SD, Spectral Volume 

• Simple, economical, and intuitive 

• Stability proofs (Jameson 2010, Vincent el al. 2011,...) 



Outline 


Review DG method 

New strong form (approximate delta functions) 
FR methods by integrating the new strong form 
Fourier and energy stability 


Conclusions 



Conservation Laws 


Conservation law 


u t + fx ~ ^ 
with initial condition 

u(x, 0) = u init (x). 

Calculate the solution ll ( X, t ) 



Legendre Polynomials 


Let P m be the space of polynomials of degree m or less. 
On I = [-1,1], for any twocontinuous functions v and w 


( v,w), 




Let the Legendre polynomial 
of degree i be denoted 
by L ( and defined by 

L ; 1 P t _ j and L ; (1) = 1. 
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Projection 


On 1= [-1,1], the projection of a function 
onto P m is 




Discretization 


1 




I 




For each cell £•, with the local coordinate g on [-1,1], 
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At time t'\ (dropping superscript n) suppose the data 
u : , are known for all j and i. 

J » 1 J 

We wish to calculate f x for ( uX +(f(uj)) x =0. 


Interface Flux 


At each interface j — 1/2, using w-_ 1/2 and w|_ 1/2 > define 
a flux fj— 1/2 (say, Roe's flux) commonfor the two adjacent cells 
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Jumps at interfaces 

On E — Ej, denote ( u,v ) E =J u(x)v(x)dx. 
Set [f] L = ft - ff and [/]* = ft - ff . 


ft = f( U j( X j- 1/2 



L 


f R =f(Uj(x j+l/2 )) 




J 


R 


Review DG Formulation 

On E, with test function <f> (degree k), 

(u h ,0) t + ((/(«a)W) = 0. 

Integrate by parts, 

(u h ,4>) t + (f<t>) dE ~ (/(«*)>&) = °- 
Allow data across cells to interact by 

iu h 4) t + (/■ I (p) dE - (f(u h ),<f> x ) - 0. 
The above is the weak form. Equivalently, 

(u h ,</>) t + f& R -fl0 L -(f(u h ),0 x ) = O. 
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Review DG Formulation 

Weak form: on E 

(Wft.fW, + /sVfi - ft 4 >l - (/(«*).&) = 0. 

With [ f] L = fl-fi and [/]* = /*'-/*, 
integrate by parts again, weobtain the strong form 

(u h ,</>) t +((f(Uh)) x ,</>) + [f]R</>R-[f] L </>L = Q- 
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The task is to eliminate (/>. 



Approximate Dirac Delta Function 

* For a fixed a on / = [ — 1,1] , let the approximate (Dirac) delta 
function to degree k at a be a linear functional on P k : 

S a ((t>) = (t>(a). 

* There exists a polynomial of degree k denoted by y ak = y a , 
i.e., y a e P k , such that 

* Proof. Set Ya=Y^i=Q b i L i- Then = 

or L m (a) = b m (L m ,L m ), or b m = L m (a)(2m + l)/2. 

That is, 
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Approximate Dirac Delta Function 




2i + 1 
2 


(- 1 )% 


and 



2i + 1 
2 
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New Strong Form 

Standard strong form 

(u h ,</>)t +((f (Uh))x,<f>) + [f]R<f>R-[f]L<f>L = Q- 
Using the approximate delta functions, 

(u h , <t>X + ((f(u h )) x , </>) + [f] R (S R , <f>) - [f] L ( S L , <p) = 0 
Using the projection onto P k , 

New strong form 

(u h ) t+ tf{(f(u h )) x )+[f] R S R -[f] L S L = 0. 



Three Members of a Family of FR Schemes 


Scheme DG 

( u h)t + <%((f ( u h))x) + 

[/]*(**.*-! + ^k) - [/kf S L , k - 1 +(-l)* ^k) = 0. 

V k ) \ L 



15 



New Strong Forms 


Strong form S 1 

(u h ) t +tf((f(u h )) x ) + [f] R S R -[f] L S L = 0. 

Strong form S2 

(u h ), + (y k {f{u h ))\ + [ f] R S R - [f] L S L = 0. 

• Derivative with no interaction : projection or interpolation; 

for form SI, interpolate via chain rule : ( f(u)) x = a(u ) u x 

• Interaction : approximate delta function, exact to degree k. 
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Energy-Stable FR (ESFR) Schemes 

Strong formS 1 and S2 for DG method (linear advection), 

0 u h ) t +a(u h ) # + Uh{ S R,k-\ + 

-m L (s L ,k - 1 + (-l)*^Z*)=0. 

ESFR schemes made simple: a, >0 

0 U h\ +a ( U h 


Key idea of the proof : Differentiate k times in % 



+ m* 




{-\) k a k 
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Reconstructing the Flux by 
Integrating the Strong Form S 1 

S 1 (u h ) t + y k {(f(u h )) x ) + [f] R S R - [f] L S L = 0. 



Flux polynomial (no interaction), i.e., 
discontinuous flux function, deg. k + 1 

Iwd (7) = ft + \\ # k kf(u h ))s)dZ 

fi p D of degree k + 1 determined by 

/ipd(“1) = /l ’ ZipdO) = /r 
an d ^-i(/ipd) = ^-1 (/ ( w /f)) 


FR: Integrate the Strong Form S 1 


SI (u h ) t + J? k ((f(u h )) x )+[fhS R -[f] L S L = 0. 


2(a). Correction function for the right boundary 


Left Radau Polynomials 



8r (£ ) = f_S R {ri)dri 
8 r =S r 

g R is of degree k + 1 : 

**(-!) = 0 , 

8r( D = l, 

i(8r) = 0. 


FR: Integrate the Strong Form S 1 


S 1 (u h ) t + y k {(f(u h )) x ) + \f] R S R - [f] L S L = 0. 


2(b). Correction function for the left boundary 


Right Radau Polynomials 




8l = ~8l 

g L is of degree k + 1 : 
8 l (~ D = l, 

,?r(D = 0, 


20 


Flux Reconstruction Form 

On E, for nonlinear conservation laws, set 

F = /iPD + \-f^L§L + [/]r8r • 

Then F is of degree k + 1 determined by 

F(-l) = f[, F(l) = f‘, 
and df_ l (F) = d?_ 1 (/( «„)). 

Also, F, =%kf(u h )) ( )+in L 8 L -Hn R 8 R . 
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Reconstructing the Flux 

Example: advection equation with k = 1. 




(a) Data (b) DG 


A Family of Fourier Stable FR Schemes 


Let g L of deg. k + 1 be defined by 

g L (-l) = l, g L { 1) = 0, 
and k additional conditions. 

For DG, 

«?-: iUl) = 0- 

For a family of stable schemes, 

«?-: 2(8l) = 0- 
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Correction functions (k = 1) 



(c) Scheme gQ a (SD) 


(d) Scheme g 2 


Second-Order FR Schemes (k = 1) 




(a) Data (b) DG 


U 



(c) Spectral Difference (SD) 


u 



(d) Scheme g 2 


Correction Functions for 
Fourier- Stable Schemes 


£dg - Rrm i 


g DG results in the DG method. 


§Ga 


k + 1 
2k + 1 


Rr.Jc + 1 + 


2k + 1 


R 


R,k 


g Ga van is hes at the k Gauss points 


§2 


2k + 1 


Rr,Jc + 1 + 


k + 1 
2k + 1 






vanishes at £ of the & + 1 Lobatto points 



Correction functions for k = 3 



-1 - 0.5 0 0.5 1 
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Fourier Analysis, k — 0 
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Correction Function 


Im 
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Fourier Analysis, k — 


_ 3£ 2 | 1 

gDG_ 4 2 4 ’ 



£__l 
2 2 ’ 


and 



n 



Spectra 
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Fourier Analysis, k — 1 


Orders of accuracy and errors 


Scheme 

Order of 
accuracy 

Coarse mesh error, 
W = 7ll 8 

Fine mesh error, 
w = ^/16 

DG 

3 

-3.2X1CT 4 -3.3xl0 -5 / 

-2.1xl0 -5 -l.lxKT 6 / 

<?Ga 

2 

-7.1x10" 4 + 2.4x10“ 3 / 

-4.6xl0“ 5 +3.1xloA 

$2 

2 

-2.5x10“ 3 +9.x10“ 3 / 

-7-lxlO” 4 +2.4xl0 -3 / 
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Fourier Analysis, k = 2 



Spectra 
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Fourier Analysis, k = 2 


Orders of accuracy and errors 


Scheme 

Order of 

Coarse mesh error, 

Fine mesh error, 


accuracy 

W = 7tl 8 

w-n! 16 

DG 

5 

-5.xl0” 7 -3.4xl0” 8 z' 

1 

X 

O 

1 

1 

FO 

• 

X 

O 

0 

#Ga 

4 

- 1. 4 xl(T 6 +8.5x10 A 

-2.2xl0” 8 +2.7xl0” 7 / 

$2 

4 

-3.2x10^ +1. 9 xlO -5 / 

1 

• 

X 

O 

1 

00 

+ 

• 

X 

0 

1 

-J 

• 
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Fourier Analysis, k = 3 


n 



Spectra 
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Re 


Fourier Analysis, k = 3 


Orders of accuracy and errors 


Scheme 

Order of 

Coarse mesh error, 

Fine mesh error, 


accuracy 

II 

£ 

OO 

II 

£ 

DG 

7 

1 

• 

X 

0 

1 

• 

X 

o 

1 

oo 

. 

-4.xlO“ 10 -2.xlO _11 j 

&Ga 

6 

-3.1xl0" 7 +1.3X10 -6 / 

-1.2xl0“ 9 +I.lxl0” 8 / 

8 2 

6 

-5.4xl0 -7 +2.3X10” 6 / 

-2.2xl0” 9 +1.9xl0“ 8 i 
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Stability 

* For solutions of degee k, if g is orthogonal t oP k _ 2 , 

then the (family) scheme is Fourier as well as energy - stable. 

* The above condition is not necessary: gLump,ch-Lo * s not 
orthogonal to any P m , but the resulting scheme is Fourier - stable. 

Open problems 

1. The collection of all g resulting in stable sc hemes remains 
to be identified . 

2. Is Fourier stability equivalent to energy stability? 
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Energy Stability 


Jameson (2010) proved that a particular SD 
scheme (recovered via FR) is energy-stable. 

Vincent, Castonguay, and Jameson (2011) proved 
energy stability for a family of FR schemes. 

Energy-stability proofs for advection and 
advection diffusion equations in ID, 2D, and 3D 
were provided by Vincent, Castonguay, Williams, 
and Jameson 

Can the current simplified proof for energy 
stability be extended to 2D, 3D, and tensor 
product cases? 



Summary 


Review DG method 

New strong forms (approximate delta functions) 

Reconstruct the flux: FR methods 

Simplified energy-stability proof 

Open problems (for grad student, 1 month of study) 

NASA Report TM-201 4-21 8 135 June 1014 (pdf) 

There is significant current research activities in FR 
methods for practical flow problems in CFD. 



Thank you for your attention. 


Questions/ Comments 
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